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We propose a mechanism which allows one to control the transmission of single electrons through a 
molecular junction. The principle utilizes the emergence of transmission sidebands when molecular 
vibrational modes are coupled to the electronic state mediating the transmission. We will show 
that if a molecule-metal junction is biased just below a molecular resonance one may induce the 
transmission of a single electron by externally exciting a vibrational mode of the molecule. The 
analysis is quite general but requires that the molecular orbital does not hybridize strongly with the 
metallic states. As an example we perform a density functional theory (DFT) analysis of a benzene 
molecule between two Au(lll) contacts and show that exciting a particular vibrational mode can 
give rise to transmission of a single electron. 

PACS numbers: 82.53.St, 34.35. +a 
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Several experiments have established that vibrational 
excitations can have a significant effect on the I- V char- 
acteristics of single-molecule junctions^— . In particular, 
the emergence of peaks in the differential conductance 
corresponding to vibrational frequencies, shows that tun- 
neling electrons interact with certain vibrational states 
of the molecule, possibly providing a means for con- 
trolling the transmission of electrons. A considerable 
amount of theoretical work have been dedicated to elu- 
cidating the effect of phonons on electronic transport in 
mesoscopic systems. Analysis of model Hamiltonians^— 
have given qualitative insight into the interaction of tun- 
neling electrons with molecular vibrations, while DFT 
based studies in conjunction with a non-equilibrium 
Green function approach show quantitative agreement 
with experiments^^—. Most efforts so far have been di- 
rected towards the influence of vibrations on transmission 
functions and I- V characteristics. In the present paper 
we will take a slightly different point of view and parti- 
tion the electronic transmission function into pieces that 
involve different vibrational excitations. We then show 
that the molecular junction may be put in a configura- 
tion where the vibrationally excited molecule allows the 
transmission of a single electron, while transmission is 
forbidden in the vibrational ground state. Controlling 
the vibrational state of the molecule, e.g. by means of a 
laser, then implies control of single electron transmission. 

The system under consideration is a molecule sand- 
wiched between two metallic leads. We assume that there 
is a single unoccupied molecular state which obtains a fi- 
nite lifetime due to hybridization with metallic states. 
This state will be referred to as the resonance and its 
position may be tuned by applying a gate voltage. If a 
bias voltage is applied to the contacts and the resonance 
is positioned in the bias window, electrons may tunnel 
through the resonance into the downstream contact and 
one will observe a current. If the molecule is weakly inter- 
acting with the metal such that the resonance is well lo- 
calized in energy, one can apply a gate voltage which sit- 



uates the resonance above the upstream chemical poten- 
tial and no current will be observed. However, if the res- 
onance couples to molecular vibrations, an off-resonant 
enhancement of transmission known as transmission side- 
bands may be observed^. In particular, if the molecule is 
initially vibrationally excited, off-resonant electrons be- 
low the resonance may tunnel though the molecule by 
absorbing a vibrational quantum of energy. This is il- 
lustrated in Fig. Q] Since transmission is only allowed 
combined with a downward vibrational transition, only 
one or a few electrons may tunnel through the contact 
and the transmission channel will be closed once the vi- 
brational mode reaches the ground state. 

Inspired by these considerations, we perform a quanti- 
tative analysis based on the model Hamiltonian&H 

H = e c\c a + ^2 M\bi + A i c K0] + b i) C 1 ) 
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where c\ is the creation operator for the lowest unoccu- 
pied molecular orbital (LUMO), c Lk and c Rk are creation 
operators for metallic states in the left and right lead 
respectively and b\ are creation operators for the vibra- 
tional normal modes of the molecule with frequencies uii. 
Thus, the electronic states of the left and right contacts 
are coupled through the molecular resonance which is 
coupled to molecular vibrations with coupling strengths 
Xi. We impose the wide band limit in which the contact 
density of states is constant in the region of the resonance 
and the resonance hybridization with metallic states is 
determined by the parameters 
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FIG. 1: A molecule between two metal contacts is represented 
by a resonant state (for example the lowest unoccupied molec- 
ular orbital) and a vibrational potential. If the molecule is ini- 
tially vibrationally excited, an electron below the resonance 
may tunnel through the molecule by absorbing a quantum of 
vibration. When the molecule is initially in its vibrational 
ground state, transmission is not possible. 



and a similar expression for the coupling to the right 
lead r^. Without the vibrational coupling, the reso- 
nance spectral function would then be a Lorentzian with 
full width at half maximum given by T = Tl + Fr- We 
will be interested in the regime T, fc^T huji, but we 
will not restrict ourselves to the classical limit r < fcgT 
where the contact current can be expressed in terms of 
rate equation o 8 ' 9 ' 14 . Instead, we consider the transmis- 
sion matrix T n (ei,Ef) for an electron with initial state 
energy ei and final state energy ef. Within scattering 
theory, the transmission matrix can be expressed in terms 
of a two-particle Green function^ which can be evaluated 
exactly in the wide band limits. The subscript n refers 
to the initial state of the oscillator and integrating out the 
final state energy, one obtains the transmission probabil- 
ity T nm (ei) that an electron with initial state energy Si 
is transmitted while the molecule makes the vibrational 
transition n — > m. 

In the appendix we have calculated the transmission 
matrix in the ground and first excited state and in Fig. 
[2] we show the transmission probability corresponding to 
four different vibrational transitions of the molecule. It 
shows that incoming electrons with energies below eQ — hu 
have a vanishing probability of transmission (Y) n Ton) 
when the molecule is in its vibrational ground state. This 
means that if a bias voltage is applied such that the Fermi 
level of the upstream contact is at ep = eo—huj no current 
will be observed when the molecule is in its vibrational 
ground state. If the molecule is vibrationally excited, 
e.g. by means of a IR laser, transmission becomes possi- 
ble through the low lying vibrational sideband (Tn and 
Tio). However, the first electron which is transmitted 
through the Tio channel induces a transition to the vi- 
brational ground state and thus closes the transmission 
channel completely. Hence, the net effect is that applying 
a short laser pulse to the molecular contact can induce 
the transmission of a few electrons. 

The distribution of the number of electrons being 
transmitted as a result of a vibrational excitation de- 
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FIG. 2: Transmission probabilities calculated from {TJ as a 
function of incoming electron energy with Fl — Tr = 0.04ftu 
and A = 0.4/kj. Below the resonance energy so the ground 
state transmission functions Too(£i) and Tbi(ei) essentially 
vanish. The insert shows the lower sideband, where stimu- 
lated emission of a vibrational quantum Tio(ei) is the domi- 
nating transmission channel. 



pends on the ratio Tn/Tio shown in Fig. [3] as a function 
of the vibrational coupling. For small coupling param- 
eters the ratio approaches zero and the first electron to 
be transmitted is therefore highly likely to induce a vi- 
brational decay and close the channel. For A <C Huj we 
thus have T\\ -C T\q which is needed for single electron 
transmission. One might worry that the inelastic trans- 
mission amplitude may become too small for anything to 
happen in this case. However, the absolute amplitude of 
sideband transmission can easily be small if the vibra- 
tional lifetime of the molecule is long. For physisorbed 
molecules the lifetime of a vibrational state is typically 
on the order of nanoseconds, whereas for example with 
Au(lll), we can use the density of states to estimate 
that ~ 40 electron hit each surface atom per picosecond 
within 0.1 eV of the Fermi level. Thus, as long as Tio is 
on the order of ~ 1 x 10 -4 there will be plenty of attempts 
to result in a single transmission event. 
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FIG. 3: Ratio of elastic and inelastic transmission probabili- 
ties in the first vibrationally excited state. At weak coupling 
the elastic transmission Tii goes to zero indicating that a 
vibrational excitation results in only a single electron being 
transmitted. 
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FIG. 4: The principle of single electron transmission. A small 
bias is applied such that the bias window just covers the res- 
onance: eVB ~ r and a gate voltage is tuned such that the 
resonance is located at: so ~ M£ ~ eVs/2 + hw. In the vibra- 
tional ground state the transmission function Too (solid line) 
is zero in the bias window. Exciting the first vibrational state 
changes the transmission function which is dominated by the 
inelastic part Tin (dashed line) in the bias window. A vibra- 
tional excitation of the molecule will thus result in a single 
electron being transmitted. 



The setup is illustrated in Fig. @] where a small bias 
voltage Vb ~ T/e has been applied and a gate voltage 
Vg has been tuned such that the position of the reso- 
nance is located Hui above the bias window. It is crucial 
that the electronic resonance has a width much smaller 
than the quantum of vibration (r <C fiw), since otherwise 
there will be a small but constant transmission probabil- 
ity when the molecule is in the vibrational ground state. 

As an example we have performed a density func- 
tional theory (DFT) study of a benzene molecule sand- 
wiched between two Au(lll) contacts. The calculations 
were performed with the code gpaw 16 ' 17 which is a real- 
space DFT code using the projector augmented wave 
method^. The contact were simulated by a three layer 
Au(lll) slab where the the top layer has been relaxed. 
We used a supercell with 12 Au atoms in each slab layer 
which were sampled by a 4 x 4 grid of K-points and 12.2 A 
of vacuum. Benzene were added with its plane parallel 
to the surface and the adsorption energy was calculated 
as a function of distance to the slab. This is shown in 
Fig. 0for PBE±2, revPBE^ . and vdW 2 ! functional. The 
PBE and revPBE functional shows weak and no bonding 
respectively whereas the vdW functional shows a 0.3 eV 
minimum at 3.6 A. The weak van der Waals bonding 
indicates that benzene is physisorbed on Au(lll) which 
means that the molecular orbitals are weakly hybridized 
with metallic states as required by vibrationally medi- 
ated transmission of single electrons. This can also be 
seen explicitly from the Kohn-Sham projected density of 
states from which we estimate the width of HOMO and 
LUMO resonances to be ~ 0.01 eV. It should be 
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FIG. 5: Adsorption energy of benzene on Au(lll) as a func- 
tion of distance to the surface, calculated with three different 
functionals. The molecular states are weakly hybridized with 
the metallic states and only the functional including the van 
der Waals interaction gives the correct adsorption well. 



noted that we assumed a LUMO state in Eq. ((T|) and 
Fig. [21 but the analysis is equally valid for transmis- 
sion of hole states mediated by the HOMO and in the 
following we will consider both types of resonance. Ben- 
zene has two degenerate HOMOs and two degenerate LU- 
MOs which have the potential to mediate transmission 
of electrons through the molecule. One of the LUMOs 
is shown in Fig. [5] and it is expected that a transient 
occupation of the orbital may induce internal forces in 
the molecule and thus couple to the vibrational modes of 
the molecule. We have performed a DFT-based normal 
mode analysis of the benzene molecule which has 36 vi- 
brational modes. There are 6 degenerate highly energetic 
modes with fkdi ~ 0.39 eV and the rest of the modes are 
evenly distributed in the interval hu>i ~ — 0.20 eV. The 
high energy modes involve the hydrogen atoms oscillat- 
ing in the plane of the molecule along the individual H-C 
bonds as shown in Fig. [6] The HOMO and LUMO states 
are expected to couple to several of the molecular vibra- 
tional states but we will focus on the high energy modes 
which have highly separated vibrational sidebands in the 
weak coupling limit. The coupling constants Ai can be 
related to the excited state potential energy surface V a 
associated with the LUMO being occupied or the HOMO 
being emptied^ 2 -: 



A, 



k dV a 



y/2 du t 



U 



(3) 



where m is the coordinate of the i'th normal mode, 
its equilibrium position in the electronic ground state, 
rrti its effective mass and V a the excited state potential 
energy surface associated with the resonanc e) 22 ' 23 . To 
obtain V a we have used the method of linear expansion 
A self-consistent field DFT 23 which allows us to calcu- 
late the excited state energies while moving the atoms 
along the high energy mode where the 6 hydrogen atoms 
move in phase (see Fig. [5]). The results for this mode 




4 




FIG. 6: Left: The lowest unoccupied orbital (LUMO) of ben- 
zene. The red and blue surfaces are positive and negative 
isosurfaces of the wavefunction. Right: The vibrational mode 
of highest energy which can be excited by a transient occu- 
pation of the LUMO. 

are \ HOMO = 27 meV and A ic/MO = 9 meV which 
should be compared with the quantum of oscillation 
huji ~ 0.39 eV. The coupling is thus rather weak and we 
obtain maximum probabilities of inelastic transmission 
of 5 x 1CT 3 and 5 x 10" 4 respectively (with T = 0.01 eV) 
at the lower vibrational sideband. However, assuming 
the vibrational lifetime to be on the order of nanosec- 
onds, the probabilities are most likely large enough that 
an elastic transmission event will occur. Furthermore, 
a small coupling constant means that the ratio Th/T±q 
becomes very small and vibrational excitation will thus 
nearly always result in only a single electron being trans- 
mitted. 

In summary, we have presented a method, which al- 
lows one to control the transmission of single electrons 
in weakly coupled molecule-metal junctions. The trans- 
mission is mediated by exciting a vibrational mode of the 
molecule while a gate voltage is tuned such that the reso- 
nant state is kept a quantum of vibrational energy above 
the bias window. It is assumed that such an excitation 
can be obtained with an external perturbation, for exam- 
ple a short laser pulse. The requirement of weak metallic 
coupling (r -c hjo) is essential since it excludes elastic 
transmission in the vibrational ground state. For small 
vibrational coupling (A <C Hco) the junction will then be 
highly reliable and always give rise to one electron being 
transmitted. For large vibrational coupling (A > hui) a 
vibrational excitation will typically result in a few elec- 
trons being transmitted due to a non-vanishing elastic 
transmission. To illustrate a quantitative approach to ob- 
tain the parameters of a real system, we have used DFT 
to calculate coupling parameters for a benzene molecule 
interacting with two gold contacts. Since benzene is 
bound by van der Waals forces and only show a weak 
hybridization with metallic states it satisfies the mini- 



mum requirement for the principle to work. However, 
there may be others reasons why this system is not well 
suited for experiments of this kind and it would be very 
interesting to investigate the principle in systems where 
transmission through single molecules with significant vi- 
brational coupling has been observed^— . 
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Appendix A 

In this appendix we will show the details of the cal- 
culations leading to Fig. [2] for a single vibrational mode 
of frequency ujq and coupling Ao . Within scattering the- 
ory the transmission matrix T n (ei,Sf) for a vibrationally 
excited state n can be expressed as£ 

T ( E - e f )-Trr n [ dTdtdS e «K««-g/fr+°/«-s"]/ft 

xG n (r,s,t), (Al) 

where 

G n (r, s, t) = 6{s)0{t) (n\c a (r - s)4 (r)c Q (t)c£ \n) (A2) 

is the two particle Green function of the vibrational state 
n. The Green function can be evaluated exactly in the 
wide band limit and the result isi£ 

G n (r,s,t) =G%(t)G^{t)e l ^(t-s) 

xe-^L n [g(f T>s , t + f^ t )}, (A3) 

where L n is the n'th Laguerre polynomial, g = 
G° R (t) = ~i9(t)e~ l{£0 ' ir/2)t/h , 

and 

f T S t = 2- e- iulot - e lul ° s + e- iUoT (l - e l " ot )(l - e lul ° s ). 

An explicit result for T n can be obtained by perfoming 
the integrals in Eq. (jAll) after a Taylor expansion of 

The result for the vibrational ground state T involves 
Lq(x) = 1 and has been calculated previously— Here we 
simply state the result which is 
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It is clear that integrating over final state energies simply produces a sum over vibrational transitions such that the 
n'th term in Eq. (|A4p represents T 0n . Calculating Ti(e.i,Ef) is a bit more involved since the integrand in (|A1|) now 
includes the first Laguerre polynomial Li{x) = 1 — x. We start by writing 



where 
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Taylor expanding the second exponential in the r integral and performing the integration gives 
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The first two terms are each others complex conjugated and the last two terms factorizcs (s and t integrals) into 
complex conjugated and the integrals can then be performed. The final result is rather complicated but consist of an 
infinite number of terms, each of which involves a delta function S(si — Ef — mfiuo), where m runs from -1 to infinity. 
We can thus obtain T 10 and Tn by collecting terms involving m = — 1 and m = respectively. The results are 
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